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ABSTRACT 

The H.E.S.S. experiment, an array of four Imaging Atmospheric Cherenkov Telescopes with high sensitivity and large field-of-view, has been 
used to search for emitters of very-high-energy (VHE, >100 GeV) y-rays along the Galactic plane, covering the region 30° < I < 60°, 280° 
< / < 330° , and -3° < b < 3°. In this continuation of the H.E.S.S. Galactic Plane Scan, a new extended VHE y-ray source was discovered 
at Q- 2 ooo=19''12'"49 ! , <5 2 ooo=+10°09'06"(HESS J1912+101). Its integral flux between l-10TeV is ~ 10% of the Crab Nebula flux in the same 
energy range. The measured energy spectrum can be described by a power law dN/dE ~ E~ r with a photon index T = 2.7 ± 0.2 S ( a t ± 0.3sys- 
HESS J 19 12+ 101 is plausibly associated with the high spin-down luminosity pulsar PSR J 1913+ 101 1. We also discuss associations with an as yet 
unconfirmed SNR candidate proposed from low frequency radio observation and/or with molecular clouds found in 13 CO data. 

Key words. Surveys - ISM: supernova remnants - ISM:individual objects:PSR J1913+1001 - Gamma rays:observations - ISM:individual ob- 
jects:HESS J1912+101 



1. Introduction 

A significant fraction of the recently discovered sources 
of very-high-energy y-rays (VHE; E > 100 GeV) within 
the Galactic plane are associated with pulsar wind 
nebulae (PWNs) of high spin-down luminosity 
sars, such as G 18.0-0.7 (lAharonian et al. 200"6dh 
MSH 15-52 (lAharonian et al. 2005bl) , or with 
panding shells of supern ova remnants (SNR s) 
as RXJ1713.7-3946 (lAharonian et al. 2007al) 



pul- 
and 
ex- 
such 
and 



RX J0852.0-4622 (lAharonian et al. 2007bl) . The observed 



Send offprint requests to: stefan.hoppe@mpi-hd.mpg.de, em- 
manuel.moulin@cea.fr 
* supported by CAPES Foundation, Ministry of Education of Brazil 
** Stanford University, HEPL & KIPAC, Stanford, CA 94305-4085, 
USA 



VHE y-ray emission from these objects provides direct ev- 
idence for the existence of particles accelerat ed to energies 
of at least ~ 10 14 eV (lAharonian et al. 2007ah . In the case 
of accelerated electrons the VHE y-rays are produced via 
inverse-Compton scattering on ambient photon fields, or, in case 
of dense environments, via non-thermal bremsstrahlung. In the 
case of accelerated hadrons, n° mesons and other hadrons are 
produced in interactions with nuclei in the interstellar medium. 
The y-rays result from the decay of these particles, in particular 
of the 7r°s. 

Most of the Galactic objects expected to emit VHE y-rays are 
related to some phase in the evolution of massive stars, either 
to their wind outflows or to PWNs and SNR shells formed 
after their collapse. Their spatial distribution is characterised by 
close clustering along t he Galactic plane (rms in latitude ~ 0.3°, 
lAharonian et al. 2006al) . A systematic survey of the Galactic 
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Fig. 1. Top, left: Image of the VHE y-ray excess from HESS J 19 1 2+ 101 , smoothed with a Gaussian profile of cr=0.13° along each 
axis. Significance contours calculated for an integration radius of $cor=0.22° are shown in black at 3.5, 4.5, 5.5 and 6.5 <x. The 
black cross indicates the best fit position of the source centroid together with its statistical errors. The dashed circle indicates the 
region from which the spectrum shown in Fig. [2] has been extracted. The positions of the pulsar PSR J1913+1011, the ROSAT 
source 1RXS J191 159.6+100814, the INTEGRAL source IGR J19140+0951 and the microquasar GRS 1915+105 are marked. The 
Galactic plane is indicated by a white dotted line. The white solid circle illustrates the 68% containment radius of the Gaussian 
brightness profile used for smoothing convolved with the point spread function of the instrument. Top, right: Image of the sur- 
face brightness (in units of MJy sr~') me asured at ~ 8.0 w ithin the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE) using the Spitzer telescope dBeniamin et al. 20031) . The white contours indicate the significance of the VHE y-ray ex- 
cess. The region marked as "B" is a complex of molecular clouds and Hll-regions. The white ellipse ill ustrates the po sition and 
extension of the SNR candidate 44.6+0.1 proposed in the Clark Lake 30.9 MHz Galactic plane survey (Kassim 1988). Bottom, 
left: Velocity profile of 13 CO(J=1->0) intensity at 110.2 GHz from the Galactic Ring Survey (GRS. lJackson et al. 20061) . integrated 
within the dashed circle shown in the top, left figure . The velocity resolution is 0.25 km s~'. The distance/velocity correspondence 
from the Galactic rotation model from Fitch et al. d 19891) is also shown. The velocity corresponding to the nominal distance of 
PSR J 19 13+ 101 1 (4.48 kpc) dManchester et al. 20051) is marked by an arrow together with the velocity range of the image shown in 
the bottom, right figure. Bottom, right: GRS 13 CO(J=l— >0) intensity in units of Kkms~'. Intensities are integrated in the velocity 
range 50 - 70 km s _1 , corresponding to a distance of ~ 5 kpc. The molecular cloud marked as 'A' is coincident with PSR J191 3+101 1. 
In region 'C no significant VHE -y-ray excess was observed so far, even though it contains molecular clouds at a similar distance. 



plane is therefore an effective approach to discovering new the High Energy Stereoscopic System (H.E.S.S.) is the VHE 
y-ray sources and source classes. Due to its large field-of-view, instrument currently best-suited to performing such a survey, 
its high sensitivity and its location in the southern hemisphere, The first phase of the survey has resulted in numerous detections 
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of VHE y-r ay emitters in the inner Galactic region \l\ < 30° 
and \b\ < 3° dAharonian et al. 2006at lAharonian et al. 2005ah . In 
2005-2007 the survey has been extended, covering the regions 
I = 280° - 330° and I = 30° - 60° . Here we present one of 
the new sources, HESS J1912+101, discovered during these 
observations, which is a good candidate for a VHE PWN such as 
HESSJ1825-137 (lAharonian et al. 2006dl). HESS J1718-385 
and HESS J1809-193 (lAharonian et al. 2007cl) . 

2. Observations and Results 

The region of interest, which incl udes the microquasar 
GRS 1915+105 dRodriguez et al. 19951). the Integral source 
IGRJ19140+0951 dHannikainen et al. 20041). the high spin - 
down luminosity pulsar PSR J 1 9 1 3 + 1 1 1 dMorris et al. 20021) . 
as well as the SNR candidate 44.6+0.1 dKassim 19881) . was 
first targeted in 2004 as part of the observational programme 
on GRS 1915+105. In 2005, dedicate d observations wer e taken 
on the supernova remnant W49B dHarris et al. 19601) : some 
of the pointings also cover the field-of-view of interest. In 
2006, the region was observed again, as part of the extended 

H. E.S.S. Galactic Plane Survey and with dedicated observa- 
tions of IGRJ19140+0951. PSRJ1913+1011 was selected for 
dedicated observations as part of a systematic study of high 
spin-down luminosity pulsars, but was not observed due to the 
planned coverage of this region in other programmes. The com- 
bined data set was investigated using the standard survey analy- 
sis (a radius of the on-source region of # cu t=0.22°, a ring back- 
ground region of radius 0.8° and hard cuts, which include a min- 
imum requirement of 200 photo electrons per sh ower im age, for 
y-ray selection) as described in Aharonian et al. (20 06al) . 

The analysis presented considers only observations pointed 
within 3° of the original source candidate apparent within the 
Galactic plane survey analysis. The mean pointing offset is 

I. 1°. After quality selection to remove data affected by unsta- 
ble weather conditions or hardware issues, the data set has an 
acceptance corrected live time (equivalent time spent at an off- 
set of 0.5°) of 20.8 h at the centre of the emission. The zenith 
angles of the observations range from 33° to 54° leading to a 
typical energy threshold of ~ 800 GeV. FigureQJtop, left) shows 
the excess count map of the 1.6° x 1.6° region around the source 
smoothed with a Gaussian profile of width 0. 13° to reduce statis- 
tical fluctuations. A clear excess of VHE y-rays is observed with 
a peak statistical significance of 7.5 cr for an integration radius of 
fl cut =0.22° . As the source was discovered as part of the H.E.S.S. 
Galactic Plane Survey, the statistical t rials associated with thi s 
survey must be taken into account (see lAharonian et al. 2006al) . 
A very conservative estimate of the number of trials which ac- 
counts for the large number of sky positions probed (~5xl0 5 ) 
leads to a post-trials significance of 5.5 cr. Fitting the uncorre- 
cted excess count map with a symmetric Gaussian profile con- 
volved with the point spread function of the instrument leads to 
a best fit position of a 2 ooo=19 /, 12'"49 J , <5 2 ooo=+10°09'6", with a 
3' statistical error in each coordinate, as indicated by the black 
cross in Fig.[U(top, left). HESS J1912+101 is clearly extended 
with an intrinsic Gaussian width of * 0.26° ± 0.03° s tat> with a fit 
probability for a Gaussian source profile of 9%. Ellipsoidal fits 
do not improve the fit probability. No further statements about 
the morphology are possible at the level of statistics available at 
the moment. 

For the spectral analysis only observations pointed within 2.0° 
of the centre of the source were used, to remove events with 
large uncertainties in reconstructed energy. The live time of 
the remaining dataset is 18.7 h. The spectrum was determined 
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Fig. 2. Differential energy spectrum of HESS J1912+101 1, ex- 
tracted from the circular region indicated in FigQ] (top, left). 
Events with energies between 0.55 and 55 TeV were used in 
the determination of the spectrum. They were binned with three 
bins per decade - the resulting number of bins is 6 - and then fit- 
ted with a power law (solid line). The two highest energy points 
of the spectrum were merged and are shown here as a 2cr upper 
limit. 

within a circular region of 0.5° radius (indicated by a dashed 
circle in FigJTJtop, left)), which represents an ~ 90% enclosure 
of the excess, chosen as a compromise between optimal sig- 
nal to noise ratio and independence of source morphology. To 
minimise spectral uncertainties, the background was estimated 
from regions with equal offset from the c entre of the field- 
of-view as described in Berge et al. d2007l) . The dataset con- 
tains a large fraction of observations targeted at other possible 
sources within the field-of-view (e.g. GRS 1915+105) and was 
taken over a time period of three years. Furthermore, due to 
the source extension, the observations of HESS J1912+101 have 
rather uneven exposure across the region and coverage over time. 
Consideration of only nearby H.E.S.S. pointing positions, use of 
the reflected-region background model and the estimation of the 
time-dependent optical response of the system, corrects for these 
variations. However, these factors lead to somewhat increased 
systematic errors with respect to typical H.E.S.S. results. Within 
the large integration circle 276 excess events were found, corre- 
sponding to a statistical significance of 7.5 cr (pre-trials) for this 
dataset used for spectral analysis. The spectrum shown in Fig. [2] 
can be described by a power law: dN/dE = Oo(£'/lTeV) _r with 
photon index T = 2.7 + 0.2 stat + 0.3 sys and flux normalisation 
O = (3.5 + 0.6 stat ± 1.0 sys ) x 10~ 12 cm" 2 s" 1 TeV" 1 . The fit 
has a ^ 2 /ndf = 6.3/4. The integral flux between 1 and 10 TeV 
is abo ut 9% of the flux of the Crab Nebula in the same energy 
range (lAharonian et al. 2006d) . 

3. Possible Associations 

The celestial region around HESS J1912+101 hosts several po- 
tential counterparts. While a superposition of two or more is 
also possible, each of the following objects could individually 
account for the observed y-ray emission: 

The most plausible counterpart candidate, both in positional 
and energetic con nection to HESS J 19 12+ 101 is the pulsar 
PSRJ1913+1011 dMorris et al. 20021) . which is slightly offset 
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from the H.E.S.S. source's best fit position by ~0.15°, spa- 
tially consistent at the ~ 3<x level. It is a rather energetic pul- 
sar with a spin-down luminosity of 2.9xl0 36 erg s _1 , a spin- 
down age of t c w 1.7 x 10 5 years, a spin period of 36 ms and 
a distance estimated fro m dispersion measurements of 4.48 kpc 
dManchester et al. 2005b . The pulsar is sufficiently energetic to 
power the H.E.S.S. source with an implied conversion effi- 
ciency from rotational energy to 1-10 TeV y-rays of e=0.5%, 
comparable to the efficiency inferred for other VHE PWN 
candidates such as HE SS J1718-385 and HESS J1809-193 
dAharonian et al. 2007c). Its distance would result in a projected 
size ofHESSJ1912+101 of ~70pc. 

These characteristics suggest an association of the VHE 
emission with the pulsar's wind n ebula, similar to othe r 
PWN associations such as VelaX (Aharonian et al. 2006bl) . 
HESSJ1825-137 dAharonian et al. 2006dl) and MSH 15-52 
dAharonian et al. 2005bl) . Compared with these examples, 
HESSJ1912+101 would be the oldest candidate for a VHE 
emitting PWN. In leptonic models of VHE y-ray production, a 
PWN is expected to also emit non-thermal X-rays. Two Chandra 
observations within the region of interest are currently avail- 
able. The first (Obs. Id 4590) was targeted on EXO 1912+097 
dLu et al. 19 96). which is coinciden t with and possibly th e 
same object as IGR J19140+0951 dHannikainen et al. 20041) . 
The overlap of this observation with HESS J1912+101 is 
marginal. The second (Obs. Id 3854) was targeted on the 
pulsar PSR J 19 13+ 101 1 and covers a larger fraction of the 
H.E.S.S. source. No X-r ay PWN was r eported at the position 
of PSRJ1913+1011 (see iGotfhelf 20041) . Existing Chandra ob- 
servations do not allow for a counterpart search on the same 
spatial scale as the H.E.S.S. source, and a detailed analy- 
sis of the Chandra data is beyond the scope of this paper. 
Despite considerably less sensitivity, the ASCA satellite with its 
larger field-of-view seems more suitable for counterpart stud- 
ies on larger scales. Two ASCA observations coincide with 
HESSJ1912+101. The first (Obs. Id 57005060) covers the 
whole H.E.S.S. source but suffers strongly from stray light con- 
tamination, the second (Obs. Id 57005070) barely covers the 
peak of the VHE emission. Therefore no ASCA data has been 
analysed for this paper. 

In the emerging picture of PWN associations the VHE emis- 
sion peak is often found offset from the pulsar position 
(e.g. HESS Jl 825- 137, MSH 15-52, HESS J1718-385 and 
HESS J1809-193). These offsets can be explained in terms 
of the proper motion of the pulsar and/or the expansion 
of the SNR/PWN into an inhomogeneous medium (see e.g. 
iBlondin et al. 2001b . In the case of HESS J1912+101 the ob- 
served offset could be explained by proper motion assuming a 
velocity of ~ 60 km s (to travel the ~ 10 pc from the centroid 
of the VHE emission to the current pulsar location in 1.7 xlO 5 
years). The hypothesis that the PWN of PSR J1913+101 1 
(and the associated SNR) is expanding into an inhomogeneous 
medium seems quite plausible as the 13 CO(J=l— »0) emission 
line measurements taken at 1 10.2 GHz indicate the existence of 
'clumpy' molecular material (in particular the area "A" as indi- 
cated in Fig.QIbottom, right)) at the pulsar position and roughly 
at its distance. Such a scenario for an asymmetric nebula wa s 
also proposed for HESS J1825-137 dAharonian et al. 2006dh . 
HESS J1912+101 and HESS J1825-137 share several similar- 
ities: 1) similar offset of the VHE emission region from the 
pulsar position; 2) similar angular size (rms size of 0.24° for 
HESS J1825-137 and 0.26° for HESS J1912+101); 3.) similar 
distances (inferred from dispersion measurements of the corre- 
sponding pulsars to be ~4kpc). The intrinsic sizes are hence 



both ~40pc. In a leptonic scenario, the cooling timescale, in 
combination with the propagation speed of the lowest energy ( 
and hence longest-lived) electrons which upscattering soft pho- 
tons to energies detectable by H.E.S.S., limit the extent of the 
emission region. As argued by Aharonian et al. (2006d) syn- 
chrotron cooling is likely to be the dominant energy-loss pro- 
cess for magnetic field strengths B > 3 juG and multi-TeV elec- 
trons. For scattering in the Thomson regime by cosmic mi- 
crowave background radiation (CMBR) photons, the parent elec- 
tron population responsible for the VHE y-rays detected by 
H.E.S.S. have energies 10-100 TeV. The magnetic field strength 
within the extende d X-ra y PWN of PSR B 1823- 13 suggested 
by Gaensler et al. d2003b is ~ 10 /jG. Considering the different 
sizes of the X-ray and VHE y-ray emission regions (5' versus 
60') for PSR B1823-13, the magnetic field strength in the outer 
VHE nebul a should be low er than in the more compact X-ray 
nebula (see Ide Jager 2008b . Here we assume similar magnetic 
field strengths for HESS J1825-137 and HESS J1912+101. The 
corresponding lifetime due to synchrotron losses of electrons up- 
scattering CMBR photons to energies E y is given by 



'synch 



(i) 



Even for a magnetic field strength of 5 yuG - a value similar to 
the interstellar magnetic field - and 1 TeV y-rays is the syn- 
chrotron cooling time scale for the upscattering electrons a fac- 
tor ~ 6 lower than the age of PSR J1913+1011 (r e = 1.7 X 10 5 
years) and stiU close to the age of PSR B1823-13 (r e = 2.1 x 10 4 
years). This provides a natural explanation for the similar sizes 
of HESS J1825-137 and HESS J1912+101 even though the 
spin-down ages of the corresponding pulsars are different by a 
factor of ~ 8. Assuming the transport process to be dominantly 
advective rather than diffusive, the necessary average flow speed 
needed to drive electrons from the position of the pulsar to the 



edge of the y-ray nebula (D 1 ™ 1 
chrotron loss time is 



40 T5k^P c ) witnin the s y n - 



D" 



v = 



'synch 



M^i^U^T^ 1 ■ (2) 



This speed is fairly large, but may not be implausible if 
the associated SNR expanded at least in part into the low- 
density phase of the interstellar medium. Assuming an 
energy-independent advection speed, equation [2] predicts an 
energy-dependent morphology for the nebula, as was observed 
for HESS J1825-137, whe re the photon index ch anges with 
distance by A T « 0.6 dAharonian et al. 2006dl) . With the 
available statistics it is not possible to extract position-resolved 
spectra for HESS J1912+101 to test this hypothesis also in this 
case. We also note that a possibly significant fraction of the 
observed VHE emission could also result from IC scattering of 
higher energy target photons, e.g. emission from the close-by 
molecular clouds and HH-regions as visible in Fig[TJ(top, right). 
This component of the VHE emission would probe lower energy 
electrons with significantly higher cooling times. Its size would 
therefore be limited not by the cooling time but by the age of 
the system. 

Whilst the spin-down powers of PSR B 1823- 13 and 
PSR J1913+1011 are very similar (~3xl0 36 erg s" 1 ), 
HESS J1825-137 is a more luminous source than 
HESS J1912+101, resulting in fairly different apparent 
VHE conversion efficiencies of 3% for the former and 0.5% 
for the latter. The spectral indices of the VHE sources 
are, however, compatible within statistical uncertainties 
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<T = 2.38 + 0.02 stat ± 0.15 sys and r = 2.7 ± 0.2 stat ± 0.3 sys ). 
Assuming a continous injection of electrons with a power law 
dN/dE w E~ a where a — 2.0 .. 2.2, a is increased to 3.0 ... 3.2 
by radiative cooling, and a VHE y-ray spectrum with a photon 
index of 2.0 to 2.1 is expected above a synchrotron break energy 
of ~ 30GeV for HESS J1912+101. In Aharonian et al. d2005cf) 
different scenarios were discussed to explain the observed softer 
photon index, including the effect of Klein-Nishina corrections 
to the IC scattering cross-section and the superposition of IC 
spectra arising from electrons injected in the past. Both these 
arguments are equally applicable to HESS J1912+101. 
Another candidate object for the observed emission is the shell 
of a hypothetical supernova remnant which would be rather 
old if associated with the birth of PSR J 19 13+ 1011 (spin-down 
age 1.7 x 10 5 years). In this scenario, the observed y-ray 
emission could be explained as the result of the interaction 
of accelerated particles from the SNR with dense molecular 
clouds. A candidate for a supernova remnant assocation is 
44.6+0.1, an yet unconfirmed SNR candidate p roposed in the 
Clark Lake 30.9 MHz Galactic plane surve y dKassim 19881) 
and further investigated by Gorham dl990l) . The estimated 
position and size of this object are illustrated by a white 
ellipse in Fig.l (top,right). 44.6+0.1 overlaps with the VHE 
emission region, but the overall match in morphology is rather 
poor. However, the resolution from the Clark Lake survey is 
insufficient for a firm conclusion about the morphology and 
only part of the remnant might be visible at 30.9 MHz. No 
SNR has been detected in X-rays. The only published keV 
X-ray sources in this region are the faint ROSAT sources 
1RXSJ191 159.6+100814, 1RXS J191254 3+103807 and 
1RXSJ191 114.8+102102 dVoges et al. 20001) . Such a lack of 
X-ray emission c ould be explained using for example the model 
of Yamazaki et al d2006l) . in which an age of the SNR of 1 .7x 10 5 
years would imply a rather low X-ray to TeV energy flux ratio 
of ~0.01 (for the energy ranges 2-10 keV and 1-10 TeV). 
In this scenario the spatial distribution of the VHE emission 
should then trace the distribution of the molecular clouds which 
act as target material. Such clouds have been observed in the 
vicinity of PSR J1913+101 1 (indicated as area "B" in Fig^ 
(top, right)) in infrared at ~ 8 fim (FigQ] (top, right)), in the 
13 CO measurements presented above (Fig Q] (bottom, right)) and 
in the GRS CS(J=2— >1) map of the region. The best distance 
estimate for G44. 3+0.1, an HII region located within this region 
and coincident with 1RXS J191 1 59.6+100814, is 3. 9 kpc (from 
HI absorption measurements, iKucharetal. 19941) . which - 
given the uncertainties of the two methods - would place it 
at a distance consistent with that of the pulsar (4.48 kpc, from 
dispersion measurements). The VHE emission peaks close to 
this dense region. There is, however, no evidence for an overall 
spatial correlation of the y-ray emission with molecular clouds 
(see e.g. area "C" in FigQ](bottom,right)). In summary, we find 
little evidence to support the scenario of an old SNR causing the 
observed VHE y-ray signal, but can not exclude this hypothesis 
on the basis of the available experimental data. 
Other celestial objects in the vicinity of HESS J1912+101 
include the microquasar GRS 1915+105 and the Integral source 
IGRJ19140+0951. GRS 1915+105 is ~l°away from the 
centroid of the VHE emi ssion. The suggested distances to 
GRS 1915+105 (1 2.5 kpc dRodriguez et al. 19951) and 6.5 kpc 
dKaiser et al. 200 5)) would correspond to a projected dis- 
tance to the centroid of HESS J1912+101 of 110 to 220pc 
and the suggested jet termination sites (IRAS 19132+1035, 
IRAS 19124+1106) are both situated outside the VHE emis- 
sion region. A possible association of GRS 1915+105 and 



HESS J 19 12+ 101 is therefore considered as unlikely. The rela- 
tive position of IGR J19140+0951, whic h is probably a Super 
Giant X-ray binary dNespoli et al. 20071) . renders a dominant 
contribution to the VHE emission from this object unlikely as 
well. 



4. Summary 

The continuation of the H.E.S.S. survey of the Galactic plane 
has resulted in the discovery of a new extended VHE y-ray 
source, HESS J1912+101. The complexity of the vicinity of 
HESSJ 1912+101 at other wavebands does not allow us to 
make a firm conclusion on the nature of the source yet. The 
most plausible scenario appears to be that the observed emis- 
sion is associated with the high spin-down luminosity pulsar 
PSRJ1913+1011, which would make HESS J1912+101 a y- 
ray PWN similar to HESS J1825-137 in both the character- 
istics of the observed y-ray emission and the energetics of the 
associated pulsar. Ultimately, an energy-dependent morphology 
change needs to be established in further H.E.S.S. observations 
to confirm the similarity with HESS J1825-137 as argued here. 
Also the PWN needs to be detected in deep X-ray observations. 
An alternative possibility is the interaction of accelerated parti- 
cles from an old SNR with molecular material in a dense region 
in the interstellar medium. Neither of the two scenarios is fully 
supported by existing multiwavelength data, but the existence of 
an energetic pulsar within the VHE emission region favours an 
association with a PWN, although the confirmation of a PWN at 
other wavebands is currently pending. 
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